Loschmidt Echo and the Local Density of States 



o 
o 

(N 

< 

in 



Oh 
I 



> 

(N 
O 

p 

od 

o 

ON 

o 



Natalia Ares and Diego A. WisniackQ 
Departamento de Fisica "J, J. Giambiagi" , FCEN, UBA, 14-28 Buenos Aires, Argentina. 

(Dated: August 5, 2009) 

Loschmidt echo (LE) is a measure of reversibility and sensitivity to perturbations of quantum 
evolutions. For weak perturbations its decay rate is given by the width of the local density of states 
(LDOS). When the perturbation is strong enough, it has been shown in chaotic systems that its 
decay is dictated by the classical Lyapunov exponent. However, several recent studies have shown 
an unexpected non-uniform decay rate as a function of the perturbation strength instead of that 
Lyapunov decay. Here we study the systematic behavior of this regime in perturbed cat maps. We 
show that some perturbations produce coherent oscillations in the width of LDOS that imprint clear 
signals of the perturbation in LE decay. We also show that if the perturbation acts in a small region 
of phase space (local perturbation) the effect is magnified and the decay is given by the width of 
the LDOS. 

PACS numbers: 05.45.Mt; 05.45.Ac; 05.45.Pq 



I. INTRODUCTION 

Reversibility and sensibility to perturbations of quan- 
tum systems are at the heart of fields of research as vast 
as quantum computation, quantum control and coherent 
transport [U, [2i S]- The interest in these subjects has 
greatly increased due to the development of experimen- 
tal techniques that enable the manipulation of a great 
number of quantum systems from photons to mesoscopic 
devices 

The suitable magnitude for measuring the stability of 
the quantum motion, as well as its irreversibility, is the 
Loschmidt Echo (LE) 0,01, defined as. 



M{t) = \{^o\expiiHt)exp{-iHot)\iPo)\^ 



(1) 



where ipo is the initial state, Hq the nonperturbed Hamil- 
tonian and H the slightly perturbed one. For conve- 
nience, we set h — 1 throughout the paper. The initial 
state ipo is usually a Gaussian wave packet and it fre- 
quently refers to an average taken over an ensemble of 
initial states randomly localized in phase space. The LE 
has been extensively studied in recent years and different 
time and per turbation strength regimes were shown to 
exist 

@,i,M[ni- As a function of time, this magnitude 

has three well known regimes. For very short times, it is 
parabolic or Gaussian, as the perturbation theory is valid 
to first order [l^]. This transient regime is followed by 
a decay that is typically exponential in chaotic systems 
Q . Finally, the LE finds a long-time saturation at values 
inversely proportional to the effective size of the Hilbert 
space. 

As a function of the strength of the perturbation, the 
decay of the LE has mainly tree different behaviors 0, [1] ■ 
When the perturbation is very small, in which a typical 
matrix element W of the perturbation is smaller than 



the mean level spacing A, the decay is Gaussian until 
M reaches its asymptotic values. If > A, this regime 
has an exponential decay, with decay rate given by the 
width cr of the local density of states (LDOS). This is 
usually called Fermi golden rule regime (FGR). Finally, 
when cr > A , with A the mean Lyapunov exponent of 
the classical system, the regime becomes independent of 
the perturbation and the decay rate is given by A. Al- 
though this Lyapunov regime seems to be universal as the 
intensive numerical studies have shown in the literature 
fiot - recent works have found a non-uniform behav- 
ior of the decay rate as a function of the perturbation 
strength. This was observed in an echo spectroscopy ex- 
periment on ultracold atoms in optical billiards and in 
theoretical studies of the kicked rotators the sawtooth 
map and Josephson flux qubits [13, 0, [3 ITgL [TtI . Sim- 
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ilar qualitative behavior was shown in Ref. [18j for local 
perturbations. In this case the authors found an oscil- 
lating regime of the decay of the LE around the value of 
the classical escape rate. 

Our main aim in this paper is to study the non-uniform 
decay rate of the LE in the strong perturbation regime 
mentioned before. We show that the origin of such behav- 
ior is related to the interplay between the width LDOS 
and the Lyapunov exponent. The LDOS is a measure of 
the action of a perturbation over a system and it gives 
the spread of the states of the unperturbed Hamiltonian 
in the bases of the perturbed one. In general, the width 
of the LDOS has a quadratic growth when the strength 
of the perturbation is small and, for increasing strength, 
the LDOS has a constant plateau a or an oscillating be- 
havior around a. This regime, which does not present a 
steady growth, is related to the finite number of states 
of the Hilbert space linked by the perturbation [To^ . We 
show that if (T ^ A we obtain the expected result that the 
decay rate of the LE is given by A when the perturbation 
is strong enough. However, if (t < A, the non-uniform 
behavior of the width of the LDOS around a seems to be 
imprinted in the decay rate of the LE, that is, the decay 
of the echo shows fluctuations that can mask the Lya- 
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punov decay. This unexpected result can be understood 
usingthe semiclassical dephasing representation of the 
LE [201 ■ We also consider local perturbations, that is, 
the case where the perturbation is applied to a window 
of the positions. These perturbations allow to change the 
width of the LDOS while the Lyapunov exponent remains 
fixed as the region of the applied perturbation is varied. 
Moreover, the use of local perturbations is motivated by 
recent semiclassical work of Ref. [l^ and by the fact 
that local perturbations are more amenable for labora- 
tory experiments. In this case we show that the smaller 
the slice of the perturbed region is, the more similar the 
decay rate of the LE and the width of the LDOS are. 

The model that we use in our study is the cat map, a 
paradigmatic system in quantum chaos that is attracting 
increasing interest in the area of quantum computation 
[2ll . [23 | . This system is completely chaotic and is very 
well suited to our work. It can be perturbed in different 
ways and the value of A can be chosen at will. Moreover, 
the perturbations can be easily applied to all (global) or 
a section (local) of the phase space. 

The paper is organized as follows. Section HIl describes 
the main characteristics of the classical and quantum cat 
maps . We introduce the perturbations used in our work 
and the LDOS, a measure of the perturbation action over 
the system. In Section IIIII we show the results of the 
decay of the LE for cat maps with different Lyapunov 
exponents using different perturbations. In this section 
we show the relation between the decay of the LE and 
the width of the LDOS. In Section HVl we make a similar 
study but for local perturbations. Finally, Section |V| is 
devoted to final remarks and conclusions. 



II. PERTURBED CAT MAPS 

Torus maps are canonical examples of classical chaotic 
systems. In this work, we will focus on cat maps, which 
are linear automorphisms of the torus that exhibit hard 
chaos. Anosov's theorem establishes that they are struc- 
turally stable, that is that orbits of a slightly perturbed 
map are conjugated to those of the unperturbed map by 
a homeomorphism. To make the paper self contained we 
present in this section the main aspects of the classical 
and quantum perturbed cat maps that are important for 
our study. 



A. Classical system 

Cat maps for a two dimensional torus considered in a 
unit periodic square are represented by matrices acting 
on the coordinates: 



V 




911 


912 




q 


p'. 




921 


922 




p. 



mod 1. (2) 



We take integer entries in the matrix G to ensure con- 
tinuity, and also Tr[G] > 2 and det[G] = 1 since the map 
is hyperbolic and conservative. 

The maximal logarithm of the eigenvalues of G defines 
its Lyapunov exponent, quantity that characterizes the 
rate of separation of infinitesimally close trajectories. We 
consider the following cat maps. 




The corresponding Lyapunov exponents are Ai = 
log(i(3 + V5)) w 0.96, As = log(2 + VS)) « 1.32 and 
A3 = log(4 -I- \/T5)) « 2.06 respectively. 

Now we introduce a perturbation of the cat map that 
is a shear in momentum that depends on coordinates 
and we take its linear part equal to zero at the origin, 
according to [23 |. 



where in particular, we considered. 



e(g) = — (cos(27rq) - cos{4:TTq)), (5) 
Zn 

being k the strength of the perturbation. We note that 
k < 0.11 for the perturbation strength to satisfy the 
Anosov theorem [2J]. We have also used a more general 
perturbation that is a shear in momentum and position 
(see below). We notice that the Lyapunov exponent of 
all the cat maps that we have considered do not change 
significantly when the mentioned perturbations are taken 
into account. 



B. Quantum system 

Quantization of perturbed cat maps has been an im- 
portant step in the quantum chaos studies because these 
simple systems allow the understandin g of many mani- 
festations of chaos in quantum systems |23l . [2^ . [25j . The 
wave function should be periodic in both position and 
momentum representation due to the periodic nature of 
the torus. Consequently, the wave functions are peri- 
odic combinations of delta functions. If in the coordi- 
nate representation the wave function has a period Aq 
with spacing Aq/N, then in the momentum representa- 
tion the period is Ap = 2TrhN/ Aq with spacing 2TTh/ Aq. 
If Ap = Ag = 1, it follows that 1 = 2nhN. Then, we 
have a Hilbert space of N dimension for a fixed value of H. 
As N takes increasing values, we reach the semiclassical 
limit. 

The propagator is obtained from the quadratic gener- 
ating function and it can be written as a N x N matrix 
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acting over a vector with components being the amph- 
tudes of each delta function. Up to a phase, it results, 



/ N 

m2 



-exp 



ittN 

512 



5229 ■ 



(6) 

The phase factor is equal to the unity if gi2 = 1. It can 
be shown that only cat maps with odd diagonal and even 
antidiagonal (or visceversa) can be quantized ^23] . In the 
case of Gi the number N of states of the Hilbert space 
must be even [2^, requirement that we have satisfied. 

In order to quantize the perturbed map, we use the 
fact that it can be written as a composite map. If P is 
the simple shear, the quantized map is thus defined by 
composing the evolution operators for the cat map and 
the shear in the same order as the classical propagator: 

U = U^U^, (7) 
being = exp[i2T:NSp{q)] and ^^f^ = e(g), so, 



k 1 

— {sin{2T:q) - -sin{4Trq)). (8) 



Now we define the scaled strength. 



C. Local Density of states 

The action of a perturbation on the eigenstate of a 
quantum system can be described by the LDOS. Let 
(f>j{k) and ipj{k) be the eigenphases and eigenfunctions 
of a perturbed cat map (fc is the strength of the per- 
turbation and j — 1,...,N). For an eigenstate i at fcoi 
considered as the unperturbed eigenstate, the LDOS is 
defined as follows, 



Afc) = ^ |(7/;,(fc)|V^,(fc„))|25(<^- [0,(fc) - 0,(fco)]) 

(13) 

where Afc = k — ko- The initial strength fco 7^ is needed 
to remove the peculiar non generic behavior of the un- 
perturbed cat map [H, . To avoid a dependence on 
the particular characteristics of the unperturbed state i 
we always make an average over this state. 

The width a of the LDOS gives us an idea of how many 
states in the base of the unperturbed system are needed 
to describe a state of the perturbed one. So, it is a good 
measure of the perturbation action over the s yste m. This 
quantity can be measured in different ways [27| . In our 
case, we take the half distance around the mean value 
that contains the 70% of the probability [13] . 



X = = kN. 

^ 2'Kh 



(9) 



Finally, we have the propagator for the perturbed cat 
map. 



N 



U{q,q) = \l -. exp 

V W12 



exp 



ittN 
L 512 



(ffiig - 2g'g + 5229' ) 



f 

— {sin{2'Kq) sm(47rg)) 

27r 2 



(10) 



To introduce more general perturbations, we consider 
shears in both momentum and positions [2^. Eq. (fTU]) 
shows that a shear in momentum is diagonal in the co- 
ordinates representation. For a shear in position, we can 
construct the propagator by changing to momentum co- 
ordinates and finding a diagonal matrix in this represen- 
tation. The change of basis matrix is the discrete Fourier 
transform F, where, 



1 f -2mlj 
=exp 
N 



N 



(11) 



In this work, we use the following shears in momentum 
and positions: 



U = U^U^F+U'^F, 
where = exp[^{cos{6np) — ism(47rp))]. 



(12) 




FIG. 1: (Color online) Width a of LDOS as a function of 
the scaled perturbation strength x for A'^ = 300 (green line), 
iV = 600 (blue line) and N = 1200 (red line). The perturba- 
tion is a shear in positions for the cat map given by G2. We 
have obtained the same results for Gi and G3. It is clear that 
a does not depend on the size of the Hilbert space and ex- 
hibits an oscillatory behavior for strong perturbations. Inset: 
a vs. X foi" the same cat map perturbed with shears in both 
momentum and positions. It remains steady for strong per- 
turbations. 

We compute the width cr for both perturbations in- 
troduced before (see Fig[T|). In the main plot of FiglU 
a is shown as a function of x when the cat map G2 is 
perturbed with a shear in momentum [Eq. [TU] . The ini- 
tial value of the strength is ko = 0.02 and the width is 
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FIG. 2: Density plots of the LDOS p(A(/),x) as a function of 
X- An average over 50 unperturbed states was done, (a) The 
perturbation is a simple shear in momentum. We can notice 
the coherence that is responsible for the oscillations found in 
the width of the LDOS. (b) The shear is in momentum and 
coordinates. In this case, no coherence is observed. 



computed for N = 300, N = 600 and N = 1200. We 
clearly see that the width of the LDOS does not depend 
on the number of states of the Hilbert space. We also no- 
tice that a has an oscillatory behavior for x > 15. In the 
small perturbation regime (x < 15), we see the quadratic 
behavior that can be obtained in a pcrturbative way [28| . 
For shears in momentum and position [Eq. (jl2p ] , we also 
compute cr as a function of % for N — 1200 [see the in- 
set of Fig. ([1])]. In contrast to the previous results, the 
width cr is constant for x <^ 15. The quadratic pcrturba- 
tive regime is present again for small perturbations. 

From the results of Fig. ([I]), it is clear that the os- 
cillatory behavior is related to the type of perturbation 
considered. In order to illustrate this behavior in more 
detail, we show the LDOS [Eq. [T3] as a function of the 
scaled strength x for the considered perturbations. We 
have made an average over 50 unperturbed states. FigO 
(a) is a density plot of p{A(j), x) when the cat map G2 is 
perturbed with a simple shear in momentum [Eq. pop ] 
and, in Fig. [2] (b), with both shears in momentum and 
positions [Eq. (HH)]. In the first case [see Fig. [2] (a)] 



and for x < 10, the distribution has a big peak with 
Lorentzian shape whose width becomes wider until its 
tails reach the limits of the Hilbert space . Then two 
small peaks appear at x ~ 10. These peaks move to the 
border and then come back to the center from both sides. 
That denotes the coherence given by this kind of pertur- 
bation. In contrast, when the shear is more general, the 
distribution is almost uniform for x > 10 and it lacks any 
coherence. We remark that the results that are shown in 
Figs. [1] and [5] for G2 are almost the same for the maps 
Gi and G3. 



III. LOSCHMIDT ECHO 

This section is devoted to the study of the decay rate 
r of the exponential decay of the LE. We compute F as a 
function of x for the cat maps Gi , G2 and G3 perturbed 
by a shear in momentum [Eq. llOj and by shears in mo- 
mentum and position [Eq. [12]. We have used N ~ 2000 
in all the computations of this Section. This value of the 
number of states of the Hilbert space allows for a good 
fit of the decay of the LE due to the fact that the long- 
time saturation occurs at a value inversely proportional 
to the effective size of the Hilbert space of the system. 
The initial states are coherent states randomly localized 
in phase space. An average over 200 initial states was 
made. 

Fig. [3] shows r as a function of x for Gi [panel (a)], 
G2 [panel (b)], and G3 [panel (c)] perturbed by a shear 
in momentum. The width a of the LDOS (solid lines) 
and the value of the Lyapunov exponents (dotted lines) 
are also plotted. In the insets we show the behavior of 
r in the weak perturbation regime. It is clear that for 
weak perturbations the decay rate of the LE and the 
width of the local density of states are equal: both of 
them follow the quadratic behavior given by the Fermi 
Golden Rule. For strong perturbations, F reaches the 
value of the Lyapunov exponent in the case of Gi , when 
the oscillations of the width of LDOS occur at greater 
value than Ai [Figl3] (a)] . However, in the last two cases, it 
exhibits an oscillatory behavior that is evidently related 
to a [FiglSfb) and (c)] . In the case of the map given by 
G2, A2 is less than the limit reached by a and we have 
just mentioned that we can observe the influence of the 
width of the local density of states in F. For G3, A3 is 
greater than that limit, and the influence seems to be 
stronger. In this case, we see in Fig. [3](c) that the decay 
F near the first peak of cr at x ~ 15 is strongly enhanced 
and the decay rate reach to F « 2A. The 2A decay was 
studied using the uniform semiclassical approach in Ref. 
p^ . This behavior deserves further investigation [29| . 

However, for the cat map given by G3, there are some 
points for 45 < x < 60 that seem to be misplaced. The 
reason for this disagreement is that the fast decrease of 
the echo caused by the large Lyapunov exponent makes 
F very difficult to determine. As we want to verify the 
behavior found for F when A is greater than the limit 
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FIG. 3; (Color online) Decay rate F of the LE (full circles) of 
a cat map perturbed with a shear in momentum as a function 
of the scaled perturbation strength x- The width a of the 
LDOS (solid lines) and the Lyapunov exponent (dotted lines) 
are also plotted, (a) The cat map is given by Gi. In this 
case Ai < (T for x > 15- Here the LDOS do not influence 
the LE decay in the strong perturbation regime, (b) The 
cat map is given by G2. In this case A2 is similar to the 
strong perturbation limit of a. Now we clearly see that the 
oscillations of the LE decay rate are strongly related to the 
ones of a. (c) The cat map is given by G3. In this case A3 > ct 
for strong perturbations. The oscillations exhibited by the LE 
decay rate follow the a ones with a larger amplitude. In the 
insets the weak perturbation regime can be clearly seen (FGR 
regime) . 

reached by a, we repeated the calculations for the cat 
map given by, 

whose Lyapunov exponent is A4 = log(8 + 3\/7)) ~ 2.77. 
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FIG. 4: (Color online) Decay rate F of the LE (full circles) of 
a cat map perturbed by a shear in momentum and position as 
a function of the scaled perturbation strength x- The width a 
of the LDOS (solid lines) and the Lyapunov exponent (dotted 
lines) are also plotted. The cat map is given by Gi (bottom), 
G2 (middle) and G3 (top). See text for details. 

The results were completely analogous. To improve the 
fitting it would be necessary to run these calculations for 
greater values of N, so the decay rate could be better 
distinguished but it would demand a really long compu- 
tational time because the computational time goes as N'^. 
In the next section we will show that local perturbations 
allow to study the case in which the Lyapunov exponent 
is much smaller that the width of the LDOS. 

In order to confirm that the oscillatory behavior ob- 
served in the decay rate of the LE is due to the infiuence 
the oscillations of the width of LDOS, we also make the 
same study as before for the perturbation [Eq. ([T^ ] that 
does not have any oscillations [see the inset of Fig. [T] 
and Fig. ^h)]. In Fig. Owe show F for d, G2 and 
G3, together with a, for this general perturbation. It is 
clear that as a has no oscillations, nor does F, and this 
behavior is the same for the different cat maps consid- 
ered. In all cases, they reach the Lyapunov value for 
strong perturbations. The last observation supports our 
conclusion about the infiuence of the oscillations of the 
width of LDOS over the decay of the LE in the Lyapunov 
regime. 

Now, we can go one step further by giving a semiclassi- 
cal interpretation of our results using the dephasing rep- 
resentation of the LE by Vanicek |20i . The LE is found 
to be, 

A^,(r,ro,i)-(p'-po)^]P (15) 

where vq is the center of the initial Gaussian wave 
packet, with dispersion A and an average momentum 
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Pq. ASs{r,ro,t) stands for the difference of action of 
the perturbed and unperturbed orbits. Now, when the 
shear caused by the perturbation is an integer number 
of 1/N (the grid of the quantum phase space) the in- 
terference between the perturbed and the unperturbed 
orbits, given by the difference of action, is constructive. 
Otherwise, it is destructive. So, when the shear is in mo- 
mentum, that condition is satisfied for certain values of 
k and consequently we observed oscillations with a wave 
length proportional to 1 /N. The wave length of the oscil- 
lations of a would be, therefore, Att'^/N in units of k [Eq. 
([5])], as we can verify in Fig[Tl where it is approximately 
47r^ in units of X- 

On the contrary, when the shear is in momentum and 
positions, the interference condition can not be satisfied 
simultaneously and the oscillations disappear. 




q 1 



IV. LOSCHMIDT ECHO FOR LOCAL 
PERTURBATIONS 

Perturbations to quantum systems are generally local 
when they come from experiments. This means that 
understanding of the general behavior of LE for local 
perturbation is an important issue. Recently, Goussev 
and co-workers have presented a semiclassical theory for 
the decay of the LE for local perturbations in billiards 
systems [l^. They show that this kind of perturba- 
tions introduce an oscillatory behavior of the decay of 
the LE around the classical escape rate. This quantity 
corresponds to the probability to escape from the billiard 
when the classical particle hits the region of the boundary 
where the perturbation is applied. 

So far, we have shown in previous sections that the 
characteristics of the perturbation modify the behavior 
of the width of the LDOS and that they shape the decay 
of the LE. With this in mind we have looked for that 
dependence in the case of local perturbations. Further- 
more, we are interested in the behavior of the decay rate 
of the LE as a function of the width of the perturbed 
region. With this aim, we have applied the already in- 
troduced shear in momentum to a window of coordinates 
of the phase space. We took the coordinate in which the 
function that describes the shear e{q) [Eq. [5] has one of 
its maxima as the center of the interval to be perturbed. 
This is shown in FigO in which we plot the scaled e{q) 
and the corresponding interval that is perturbed. 

Let us first see the influence of a local perturbation in 
the behavior of LDOS. In Figl6]we exhibit the distribu- 
tion p{A(j), x) when the perturbation is applied in a region 
of w = 0.4. It is clear that the distribution is bounded 
and we can see a coherent behavior similar to the one 
observed in Fig[2](a). We notice that p(A(/), x) does not 
present the same symmetry as around A(j) = as is ob- 
served in FigUl (a) when the perturbation is global. This 
is related with the position go as center of the perturbed 
region 

In Figs. [7] and Fig[5]we show the behavior of the decay 



FIG. 5; Schematic figure showing the local perturbation that 
was used. The scaled shear [27re(q)/fc] {k is the strength of 
the perturbation) is plotted as a function of the q coordinate. 
qo indicates the center of the interval that was perturbed and 
w its width. 
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FIG. 6: Density plot of the LDOS p{A.(j>, x) as a function of 
the strength x when the perturbation is a simple shear in 
momentum applied in a region of width w = 0.4 (see Fig. 5). 
It can be clearly seen that it is restricted to the perturbed 
phases and that it has a coherent behavior. 



of the LE for a local perturbation. All the computations 
in this section were done with N ~ 800 (the total size 
of the Hilbert space) . We have considered several widths 
where the perturbation was applied. As an example, in 
Figl7]we show the decay F of the LE with symbols for: y 
[w = 0.4), A (w = 0.3), + (w = 0.2) and * (w = 0.1). We 
notice that, although the oscillatory behavior is present 
for all the values oiw considered, its amplitude is smaller 
when w decreases. Moreover, the value around the oscil- 
lation takes place also becomes smaller as w becomes 
thinner. To conclude the analysis, we also compute the 
width (T of the LDOS for all the different apertures con- 
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FIG. 7: (Color online) Decay F of the LE of a cat map per- 
turbed locally with a shear in momentum as a function of the 
scaled perturbation strength x- The symbols correspond to 
the following width of the perturbed region: \/ {w — 0.4), A 
{w = 0.3), + {w ^ 0.2) and * {w = 0.1). The width a of the 
LDOS is also plotted with solid lines. It is noticeable that, 
the smaller is w, the more similar to cr F is. In all the cases 
considered, the oscillations are equal for both magnitudes. 



sidered. The width a is shown in Fig[7] with solid lines. 
It is clear that for w — 0.1 and w — 0.2, the behavior of F 
and a are identical. That is, the decay of the LE is given 
by the width of the LDOS. This result is in agreement 
with the semiclassical study for billiard systems of Ref. 
(TH - For bigger w {w = 0.3 and w = 0.4 in our case), the 
oscillations are alike, but its amplitude is not so similar. 
So, in the case of a local perturbation, the influence of 
the LDOS width over F is enhanced. 

In FiglS] we show the mean value of the decay of the 
LE F after the quadratic growth [x > 20, see Fig. [7] 
as a function of the width w of the perturbation. We 
can clearly see that f w 2w for w < 0.4. We remark 
that our result agrees with the one obtained for billiard 
systems in Ref. [ISj because the width w corresponds to 
the classical escape rate in our system. 



V. FINAL REMARKS 

In this work we have studied the decay of the LE 
for strong perturbations, that is, the regime after the 
quadratic growth of the decay (FOR regime). We have 
shown that the expected monotonic regime in which the 
decay is given by the Lyapunov exponent of the classical 
system is not universal. In fact, the behavior of the LE in 
this regime depends on the ratio between the mean value 
of the width of LDOS a and the value of the Lyapunov 
exponent. 

The general idea is the following. When a system is 
perturbed and the perturbation affects a finite number of 



FIG. 8: Mean value F of decay rate of the LE as a function of 
the width w of the applied local perturbation. The function 
2w is also plotted with dotted line . See text for details. 



states of the Hilbert space, instead of growing monoton- 
ically, the width of the LDOS reaches a plateau a or an 
non-uniform regime around a. If a/X » 1, the decay 
rate of the LE is given by A and does not depend on the 
perturbation. But, if a/X < 1 and the width a of the 
LDOS has a non uniform behavior, a non-monotonicity 
is imprinted in the decay of the LE. This regime was re- 
cently reported in the literature [H, [ij, [TgI. llTj but not 
the conditions for it to happen. 

We have illustrated the previous ideas using cat maps 
with different Lyapunov exponents and perturbations. 
When the cat map is perturbed with a simple shear in 
momentum, we have found pronounced oscillations of the 
width of the LDOS, that are a consequence of the coher- 
ence of perturbed and unperturbed orbits that this kind 
of perturbation introduces. Significantly, we have shown 
that this oscillating regime was also present in the decay 
rate of the LE when the value around the LDOS width 
oscillates is close enough to the Lyapunov exponent of 
the classical system. Moreover, when the perturbation is 
more general, the coherence between orbits is destroyed, 
so the width of the LDOS reaches a plateau and so does 
the LE decay. 

We have also studied the decay of the LE in cat maps 
when the perturbation is local [l8| . The important point 
in this case is that varying the width of the region in 
which the perturbation is applied we can change the 
width of the LDOS while the Lyapunov exponent re- 
mains fixed. With this in mind, we have perturbed a 
range of coordinates with a simple shear in momentum. 
We have found that the decay rate is almost identical to 
the width of the LDOS when the range of the perturbed 
coordinates is very small. Moreover, we have shown that 
the amplitude of the oscillations of the LDOS decreases 
with the width of the applied perturbation and the mean 
value is given by the classical decay rate. 
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